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DYNAMIC  RESPONSE  OF  AN  INSONIFIED  SONAR  WINDOW  INTERACTING 
WITH  A  TONPILZ  TRANSDUCER  ARRAY 


1.  INTRODUCTION 


One  very  typical  sonar  design  consists  of  a  soft  material  that  is  in  contact  with  an  array  of 
Tonpilz  transducers.  This  soft  material  is  usually  called  a  window  (or  sometimes  a  screen).  The 
window  serves  several  purposes:  it  protects  the  transducers  from  water,  impact,  and  debris;  it 
attenuates  nonacoustic  energy;  it  provides  a  safe  covering  to  the  sonar  system  when  it  is  shipped 
or  handled;  and  it  helps  to  minimize  turbulence  and  hydrodynamic  drag.  Tonpilz  transducers  are 
a  reliable  transducer  design  that  has  been  refined  for  many  years.  Single  resonant  Tonpilz 
transducers  typically  consist  of  a  head  and  tail  mass  separated  by  piezoelectric  stack  that  emits  a 
voltage  when  it  is  subjected  to  an  applied  force.  Tonpilz  transducers  are  a  useful  design  because 
they  can  operate  in  transmit  (active)  and  receive  (passive)  mode.  They  are  functional  in 
environments  where  higher  target  signal-to-noise  ratio  is  needed  than  can  be  obtained  from 
passive  sonar  only. 

Tonpilz  transducers  have  been  studied  in  literature  for  many  years.  Basic  design  guidelines 
exist  in  textbooks.1  General  transducer  modeling  techniques  have  been  previously  developed.2 
Finite  element  models  of  the  transducers  yield  the  mechanical  behavior  when  it  is  in  contact  with 
a  heavy  fluid  field3  and  in  contact  with  air.4  Equivalent  circuit  models  have  been  derived  and 
analyzed.5,6,7  Sonar  window  models  have  been  previously  studied.  Insertion  loss  and  echo 
reduction  measurements  are  common  in  the  literature.8,9,10  Theoretical  analysis  has  been 
developed  based  on  thick  plate  theory.9,11,12  Structurally  stiffening  the  window  for  hull 
applications  has  also  been  investigated.13  References  1  through  7  are  transducer  modeling,  and 
references  8  through  1 3  are  window  modeling  and  testing.  It  is  noted  that  there  is  no  known 
analytical  model  that  couples  the  Tonpilz  transducers  to  a  fully  elastic  sonar  window. 

This  report  derives  and  evaluates  an  analytical  model  of  a  fully  elastic  sonar  window  in 
contact  with  an  array  of  Tonpilz  transducers.  The  window  is  insonified  by  a  planewave  at 
varying  arrival  angles  and  frequencies.  This  model  is  intended  for  broadband  frequency  analysis 
of  a  sonar  system  when  there  is  significant  interaction  between  the  window  and  the  array  of 
transducers.  The  formulation  of  the  problem  begins  with  elasticity  theory,  that  models  the 
motion  in  the  window  as  a  combination  of  dilatational  and  shear  waves.  From  this  theory, 
expressions  for  plate  displacements  in  the  normal  and  tangential  direction  are  obtained.  The 
displacements  are  then  inserted  into  stress  relationships  and  these  equations  are  set  equal  to  the 
forces  acting  on  the  structure  by  the  transducers  and  the  pressure  of  the  incoming  acoustic  wave. 
The  problem  is  then  written  as  an  algebraic  system  of  equations,  in  matrix  form,  where  the  left- 
hand  terms  represent  the  zero-order  window  dynamics  and  are  equal  to  an  infinite  number  of 


right-hand  terms  that  represent  the  forces  acting  on  the  structure.  Rewriting  this  zero-order 
dynamic  term,  by  increasing  and  decreasing  the  index,  results  in  an  expression  for  the  higher- 
order  modes  interacting  with  the  applied  forces.  The  integer  shift  property  is  then  applied  to  the 
right-hand  side  of  all  of  the  terms,  resulting  in  an  infinite  set  of  equations  that  model  the  wave 
propagation  coefficients  of  all  the  modes  of  the  structure.  This  set  of  equations  is  truncated  to  a 
finite  number  of  terms,  and  solutions  to  the  displacement  fields  are  calculated.  The  transducer 
output  is  written  as  a  function  of  the  displacement  field  at  the  bottom  of  the  window,  and  this 
term  is  calculated  as  a  transfer  function  of  voltage  divided  by  applied  pressure  versus  arrival 
angle  and  frequency.  A  numerical  example  is  included  where  the  array  beam  pattern  is 
generated  at  several  frequencies  and  the  results  are  discussed. 
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2.  SYSTEM  MODEL 


The  system  model  is  that  of  a  sonar  window  attached  to  an  array  of  Tonpilz  transducers,  as 
shown  in  figure  1.  This  mechanical  problem  is  analytically  modeled  by  assuming  the  sonar 
window  is  a  fully  elastic  plate  and  the  Tonpilz  transducers  are  discrete  mass-spring-mass 
systems,  as  shown  in  figure  2.  The  plate  (or  sonar  window)  has  a  thickness  of  h  (m)  and  is 
loaded  on  the  top  surface  with  a  normal  (pressure)  forcing  function.  The  transducers  on  the 
bottom  of  the  window  are  equally  spaced  at  a  distance  of  L  (m)  in  the  jc-direction  and  each  has  a 
head  mass  per  unit  length  Mh  (kg/m),  tail  mass  per  unit  length  Mr  (kg/m),  and  stiffness  per  unit 
length  Ks  (N/m2).  The  model  uses  the  following  assumptions:  (1)  the  forcing  function  acting  on 
the  plate  is  a  planewave  at  a  definite  wavenumber  and  frequency;  (2)  motion  is  normal  and 
tangential  to  the  plate  in  one  direction  (two-dimensional  system);  (3)  the  plate  has  infinite  spatial 
extent  in  the  ^-direction;  (4)  the  head  mass  has  translational  degrees  of  freedom  in  the  x-  and 
z-directions;  (5)  the  tail  mass  and  the  spring  have  translational  degrees  of  freedom  in  the 
z-direction;  (6)  the  particle  motion  is  linear;  and  (7)  the  fluid  medium  is  lossless. 


Tonpilz  Transducers 


Figure  1.  Sonar  Window  with  Tonpilz  Transducers 


3 


Figure  2.  Model  of  Sonar  Window  with  Tonpilz  Transducers 


The  motion  of  the  elastic  plate  is  governed  by  the  equation 


pV2u(x,y,z,t)+{A  +  /v)VV  •  u  (x,y,z,t)=p 


d~u(.x,y,z,/) 
dt2 


(i) 


where  p  is  the  density  (kg/m3),  A  and  p  are  the  complex  Lame  constants  (N/nr ),  /  is  time 
(seconds),  •  denotes  a  vector  dot  product,  and  u(jc ,y,z,t)  is  the  three-dimensional  Cartesian 
coordinate  displacement  vector  and  is  written  as 


ux(x,y,z,t) 

y/x(x,y,z,t )' 

u  {x,y,z,t)  =  < 

uy(x,y,z,t) 

>  =  V^(jc,y,z,t)+V  x< 

y/y(x,y,z,t ) 

uz(x,y,z,t) 

yz(x,y,z,t ) 
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(2) 


where  <f>  is  a  dilatational  scalar  potential,  V  is  the  gradient  operator,  x  denotes  a  vector 
crossproduct,  and  is  an  equivoluminal  vector  potential.  The  formulation  is  now  condensed 
into  a  two-dimensional  problem;  thus,  uy=  0  and  d(-)/  dy  =  0  .  Expanding  equation  (2)  and 
breaking  the  displacement  vector  into  its  individual  nonzero  terms  yields 


ux(x,z,t) 


and 


dj(x,z,  t) 

dx  dz 


u:(x,z,t) 


dj(x,z,  t) 

dz  dx 


(3) 


(4) 


Equations  (3)  and  (4)  are  next  inserted  into  equation  (1),  which  results  in  two  decoupled  wave 
equations  given  by 


CjV20(x,z,t) 


2 

d  </>(x,z,t) 
~d?~ 


(5) 


and 


cjv2i//  (x,z,t)  = 


d2y/v{x,z,t ) 
~d? 


(6) 


where  equation  (5)  corresponds  to  the  dilatational  component  and  equation  (6)  corresponds  to  the 
shear  component  of  the  displacement  field.  Correspondingly,  the  constants  cj  and  cs  are  the 
complex  dilatational  and  shear  wavespeeds  (m/s),  respectively,  and  are  determined  by 


cd 


A  +  2  fu 
P 


(7) 


and 


c 


s 


(8) 
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The  equations  of  motion  are  formulated  as  a  boundary  value  problem  using  four  equations  of 
stress  written  in  terms  of  the  plates’  displacements  and  corresponding  forcing  functions.  The 
plate  is  loaded  by  a  normal  (pressure)  forcing  function  as  shown  in  figure  2;  thus,  the  normal 
stress  at  z  =  b  is  written  using  a  force  balance  between  the  pressure  in  the  fluid  and  the  plate  as 


”zz 


,  ,  .  ...  -  duz(x,b,t)  duK(x,b,t) 

( x,b,t )  =  (A  +  2 p) — - -  +  A — - =  -pa(b,x,t) , 


dz 


dx 


(9) 


where  pa(x,t)  is  the  pressure  field  in  contact  with  the  top  of  the  plate  (N/m").  The  tangential 
stress  on  the  top  of  the  plate  is  modeled  as  a  free  boundary  condition  and  is  written  as 


dux(x,b,t)  duz(x,b,t ) 


dz 


dx 


=  0. 


(10) 


The  plate  is  loaded  by  the  forces  in  the  Tonpilz  transducers’  head  masses  acting  on  the  bottom  of 
the  plate,  thus,  the  normal  stress  at  z  =  a  is 


rzz(x,a,t)  =  (A  +  2p) 


duz(x,a,t)  5  dux{x,a,t )  "  00 


dz 


■  +  A- 


dx 


=  ^fz(a^)d(x-nL), 


Yl  —  — QO 


(11) 


where  /.(«,/)  is  the  force  per  unit  length  that  each  Tonpilz  transducer  exerts  on  the  plate  in  the 
z-direction  and  d{x-nL)  is  the  Dirac  delta  function  that  distributes  the  transducer  forces  discretely 
and  periodically.  Similarly,  the  tangential  stress  on  the  bottom  of  the  plate  is 


rzx(x,a,t)  =  p 


dux(x,a,t)  duz(x,a,t) 


dz 


dx 


n- oo 


=  Y,fx(ad)d(x-nL), 


j  n=-oo 


(12) 


where  fx(a,t)  is  the  force  per  unit  length  that  each  Tonpilz  transducer  exerts  on  the  plate  in  the 
x-direction. 

The  acoustic  pressure  in  the  fluid  medium  is  governed  by  the  two-dimensional  wave 
equation  and  is  written  in  Cartesian  coordinates  as 


r?2/7a(x,z,/)  d2pa(x,z,t)  1  ^2pq(x,z,Q 
dz2  dx2  c2  dt2 


(13) 
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where  pa(x,z,t)  is  the  pressure  (N/m2)  and  Cf  is  the  real  valued  compressional  wavespeed  in 

the  fluid  (m/s).  The  interface  between  the  fluid  and  solid  surface  at  z  =  b  satisfies  the  linear 
momentum  equation,  which  relates  the  acceleration  of  the  plate  surface  to  the  spatial  gradient  of 
the  pressure  field  by 


Pf 


2 

d  uz(x,b,t) 

7? 


dz 


(14) 


where  pf  is  the  density  of  the  fluid  (kg/m3). 

The  system  output  is  a  transfer  function  of  transducer  voltage  divided  by  incident  pressure, 
commonly  called  receive  voltage  sensitivity  (rvs).  This  is  equal  to 


rvs(x,t)  =  g33  t  fs(x,t). 


(15) 


where  g33  is  the  material  constant  that  relates  the  mechanical  force  to  electrical  output 
(volts  m/N),  t  is  the  thickness  of  each  piezoelectric  piece  from  the  transducer  stack,  and  fs(x,t ) 
is  the  force  in  the  transducer,  which,  in  the  case  of  this  model,  is  the  force  across  the  spring. 
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3.  TRANSFORMATION  INTO  THE  WAVENUMBER-FREQUENCY  DOMAIN 


The  problem  is  now  transformed  into  the  wavenumber-frequency  domain  (kco)  by  using  the 
functional  form,  where  the  field  variables  are  equal  to  an  unknown  function  in  the  z-direction 
times  an  exponential  function  in  the  x-direction  times  an  exponential  function  in  time.  The 
function  in  the  z-direction  is  also  a  function  of  wavenumber  and  frequency.  The  displacements 
become 

ux(x,z,  t)=  Ux  ( k,z ,  ru)exp(i&x)exp(i  cot)  (16) 

and 

uz(x,z,t )  =Uz(k,z,o))exp{\kx)Qxp(\cot) ,  (17) 

and  the  pressure  is  written  as 

pa(x,z,t )  =  Pa(z,k,co)  exp(i kx)  exp(i«/) ,  (18) 


where  i  =  V-T,  k  is  the  wavenumber  with  respect  to  the  x-axis  (rad/m)  and  co  is  frequency 
(rad/s). 

Inserting  equations  (16)  and  (17)  into  equations  (2)  through  (6)  and  solving  the  differential 
equations  gives  the  transformed  displacement  field  in  the  x-direction  as 


Ux  ( k ,  z,  a>)=  A{k,(o)\k  exp(i  az)  +  B(k,co)\k  exp(-iorz)  + 

-  C(k,  co)\j3  exp(i  J3z)  +  D(k,co)\j3  exp(-i/?z)  ( 1 9) 

and  the  transformed  displacement  field  in  the  z-direction  as 

Uz(k,z,o))=  A(k,co)'iaexp(\az)-B(k,{o)iaexp(-\az)  + 

C(k,  a))\k  exp(i  J3z)+D(k,  co)\k  exp(-i  fiz) ,  (20) 

where  A(k,co),  B(k,  co),  C(k,co),  and  D{k,(o)  are  unknown  complex  wave  propagation  coefficients 
of  the  plate,  a  is  the  modified  wavenumber  (rad/m)  associated  with  the  dilatational  wave  and  is 
expressed  as 


a  = 


- k 2 


(21) 
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where  kj  is  the  dilatational  wavenumber  and  is  equal  to  col cd,p  is  the  modified  wavenumber 
(rad/m)  associated  with  the  shear  wave  and  is  expressed  as 


(22) 


where  ks  is  the  shear  wavenumber  (rad/m)  and  is  equal  to  co/  cs . 

Inserting  equation  (18)  into  (13)  and  solving  the  resulting  ordinary  differential  equation  yields 
the  transformed  pressure  field  as 


Pa(z,k,(o)  =  M(k,(o)exp(\yz)  +  P,((o)exp(-\yz), 


(23) 


where  the  first  term  on  the  right-hand  side  represents  the  outgoing  reradiated  pressure  field 
caused  by  the  plate  displacement  and  the  second  term  represents  the  incoming  applied  incident 
pressure  field  (the  forcing  function)  acting  on  the  structure.  In  equation  (23),  y  is  the  modified 
wavenumber  (rad/m)  associated  with  the  fluid  and  is  expressed  as 


(24) 


where  yis  purely  real  or  imaginary,  depending  on  the  sign  of  the  argument  under  the  radical. 
When  the  sign  of  the  argument  is  positive,  the  analysis  is  in  the  acoustic  cone;  when  the  sign  of 
the  argument  is  negative,  the  analysis  is  in  the  nonacoustic  region.  For  acoustic  sonar  response, 
the  analysis  is  typically  studied  in  the  acoustic  cone.  The  relationship  between  the  arrival  angle 
of  an  acoustic  wave  and  its  wavenumber  is 


(25) 


where  #is  the  arrival  angle  of  an  incoming  acoustic  wave  (rad)  with  0  corresponding  to 
broadside  excitation. 

The  forces  exerted  by  the  Tonpilz  transducers  can  be  determined  with  a  dynamical  model 
of  a  mass-spring-mass  system.  The  force  per  unit  length  in  the  z-direction  for  each  transducer  is 
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/*(«.*)  = 


U:  ( k,a ,  c o )  exp(ifcx)  exp(i<y/) 


co  M M y  —  co~  K §  ( M j  +  M yi  ) 

K  $  —  co  ~  M  y 

=  F:(co)U  ,(k,a,co)e\p(ikx)e\p(\cot),  (26) 

and  the  force  per  unit  length  in  the  ^-direction  for  each  transducer  is 


fx(a,t)  =  -co1  M H  Ux(k,a,co)e\p(ikx)exp(\cot) 

=  Fx(co)Ux(k,a,co)exp(ikx)exp(icot),  (27) 

as  the  transducer  spring  constant  is  zero  in  the  horizontal  direction.  It  is  noted  that  if  the 
transducer  is  a  double  or  triple  resonant  type,  the  expressions  given  in  equations  (26)  and  (27) 
can  be  changed  to  reflect  these  dynamic  effects. 

Finally,  the  output  of  the  transducer  can  also  be  transferred  into  wavenumber-frequency 
space  as 


rvs(jc,/)  =  /?FS(&,<y)exp(ifar)exp(i&>/) .  (28) 

Substituting  the  spring  dynamics  yields  the  receive  voltage  sensitivity  in  terms  of  the 
displacement  at  the  bottom  of  the  sonar  window  as 


RVS(k,co) 


co2 

2 

K$  —  co  Mj 


g33  t  (Ks  /w)Uz(k,a,co), 


(29) 


where  w  is  the  width  of  the  transducer  head  in  the  x-direction  (m). 

The  four  boundary  value  equations  (equations  (9)  through  (12))  are  now  rewritten  in  the 
wavenumber-frequency  domain  using  equations  (16)  through  (27).  They  become 


Tzz(k,b,co)  =  (A  +  2/j) 


dUz(k,b,co)  1  dUx(k,b,co) 


dz 


■  +  A- 


dx 


2  ' 
03  Pf 

i  r 


Uz(k,b,co)  =  -2PI{co ),  (30) 
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Tzx(k,b,a>)  -  /u 


dU  x(k,b,co)  dU  z(k,b,co) 


dz 


■  +  ■ 


=  0, 


(31) 


7 lz(k,a,a>)  =  (A  +  2ju)  +  x  =  F2(<y)  '^U z{k,a,co)8{x-nL) ,  (32) 


dz 


dx 


and 


Tzx(k,a,co)  -  /j 


dUx{k,a,co)  dU,(k,a,a)) 


dz 


dx 


n=oo 


-  Fx{(d)  x(k,a,co)S{x -  nL).  (33) 


n=-oo 
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4.  ANALYTICAL  SOLUTION 


The  analytical  solution  is  now  calculated.  This  begins  by  transforming  the  right-hand  side 
of  equations  (32)  and  (33)  so  that  the  delta  functions  are  absorbed  into  the  displacement 
functions.  This  transformation  is  derived  and  discussed  elsewhere,14  and  produces  the  following 
equations 


n~°°  F  (ri\  17=00  1  nr) 

Fz{co)  YiUz(k,a,a>)S(x-nL)  =  -*f-!-  X  Uz(k  +  —  ,a,a>),  (34) 

n=- <x>  ^  n=-cc  ^ 


and  similarly 


Fx{co)  "^Ux{k,a,co)5{x-nL)  =  ^f^  X  Ux(k  +  ^,a,a>).  (35) 

n=-oo  n=- oo 


Next  the  functional  form  of  the  displacements  from  equations  (19)  and  (20)  are  inserted 
into  equations  (30)  through  (35)  and  the  following  algebraic  matrix  equation  is  obtained 


[  A <0)  (A" )]  {y (0)  ( Ar)}= 


X  [F(”>(/r  +  —^—)]{y^n)(k 


+ 


2nn 

~L 


)}  +  P. 


(36) 


where  [A^°*(£)]  is  a  four  by  four  matrix  that  models  the  dynamics  of  the  plate  for  n  =  0, 
{y(0)(£)}  is  the  four  by  one  vector  of  wave  propagation  coefficients  for  n  =  0, 

[F(w)(^  +  Inn  /  L)]  is  the  four  by  four  matrix  that  represents  the  periodic  transducer  loading  on 
the  structure  for  «th  mode,  {y  (n)(k  +  Inn  I L)}  is  the  four  by  one  vector  of  wave  propagation 
coefficients  for  nth  mode,  and  p  is  the  four  by  one  vector  that  models  the  plane  wave  excitation. 
The  entries  of  the  matrices  and  vectors  in  equation  (36)  are  listed  in  the  appendix.  To  facilitate  a 
solution  to  the  problem,  index  shifting  is  employed.  The  integer  shift  property  of  an  infinite 
summation  is  applied  to  equation  (36),  which,  because  of  the  summation  from  minus  infinity  to 
positive  infinity,  results  in 
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Lt  L/ 

2  n{n  +  m) 


2n(n  +  m) 


=  Y  [F^+m\k  +  ^^^)]{y(n+m\k  +  ^^  '"/)}  + 


n=- oo 
m- oo 


=  X  [F(">(*  +  ^)]jy'”>(*  + 


2;rm, 


(ot)/ 


2;ztw 


L 

[p  m-  0 


[  p  m  =  0 
0  m  *  0 


(37) 


Z,  I  0  m  *  0. 


where  the  0-term  is  a  four  by  one  vector  whose  entries  are  zeros.  Once  the  [A*"'*]  matrix  is 
integer-indexed  and  the  transducer  load  matrix  indexes  have  been  shifted,  the  system  equations 
can  be  rewritten  using  all  the  n-indexed  modes  as 


Ay  =  Fy  +  P  , 


(38) 


where  A  is  a  block-diagonal  matrix  and  is  equal  to 


A  = 


[A(-D  0 

Lj 

o  [A(0)(*)] 

0  0 


0 


0 

[A  W(k  +  ^)] 

Lj 


(39) 


F  is  a  rank  deficient,  block-partitioned  matrix  and  is  written  as 


F  = 


[F(-')(k-^)]  [F(0)(*)]  [F  (X\k  +  ~)\ 
[F  [F(0)(A:)]  [F  (X\k  +  ~-)\ 

[F(-1)(A:  [F(0)(*)]  [F{l\k  +  lj-)] 


(40) 
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P  is  the  planewave  load  vector 


P  = 


0T 


(41) 


and  y  is  the  wave  propagation  coefficient  vector  that  contains  all  the  unknown  indexed 
coefficients  as 


<yH)(*-f»T  (y<0,(*))T  {y">(*+f)fT 


(42) 


where  the  unknown  wave  propagation  coefficients  are  contained  in  the  equations  as 


{y(0)(*)}  =  {zl(*,<y)  B(k,co)  C(k,<o)  D{k,co)}J 
=  [4(0)  fi(0)  C(0)  D(0)  }T  . 


(43) 


The  0-term  in  equation  (39)  is  a  four  by  four  matrix  whose  entries  are  all  zeros  and  the  0-term  in 
equation  (41 )  is  a  four  by  one  vector  whose  entries  are  all  zeros.  Equation  (38)  is  assembled,  and 
the  wave-propagation  coefficients  that  reside  in  the  y  vector  are  found  by 

y  =  [A  -F]-1  P .  (44) 

When  the  coefficients  are  determined,  the  displacements  of  the  system  can  be  calculated  using 
equations  ( 1 9)  and  (20)  and  the  n  =  0  wave  propagation  coefficients. 
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5.  MODEL  VALIDATION 


The  sonar  window  —  Tonpilz  transducer  interaction  model  —  can  be  compared  and 
validated  for  a  thin  plate  at  low  frequencies  using  a  Bemoulli-Euler  model  that  has  been 
previously  developed.14,15,16  The  stiffeners  in  these  models  are  replaced  by  Tonpilz  transducer 
dynamics  so  that  the  validation  example  corresponds  to  the  thick  plate  model  developed  in 
sections  2  through  4.  The  thin  plate  model  has  one  degree  of  freedom  that  is  the  displacement  in 
the  z-direction.  This  equation  is  written  as 


UAkM 

P,{co) 


2  T(k,co) 


(&>)  _ . ,  .  F.(co)  Inn  x 

+  ^f-LT(k,co) - £  T{k  +  —  ,co) 


(45) 


where 


T(k,co) 


Dk 4 


-1 


2 

-  phoj  + 


2  V 

PfCO 


i  r 


(46) 


and 


D  = 


B )3 

12(1  -l>2)  ’ 


(47) 


where  E  is  Young’s  modulus  (N/m2)  and  u  is  Poisson’s  ratio  (dimensionless).  Once  the  normal 
displacement  is  known,  the  transducer  output  is  determined  using  equation  (29). 

Figure  3  is  a  plot  of  the  transfer  function  of  transducer  voltage  divided  by  input  pressure 
versus  wavenumber  at  a  frequency  of  1 00  Hz.  This  extremely  low  frequency  was  chosen 
because  it  is  a  value  at  which  the  two  models  should  theoretically  agree.  Additionally,  the 
comparison  is  made  versus  wavenumber  rather  than  arrival  angle  so  that  the  higher  wavenumber 
dynamics  are  included  in  the  comparison.  This  example  was  generated  with  the  following 
system  parameters:  window  thickness  h  is  0.005  m,  window  density  p  is  1200  kg/m3.  Lame 
constant  A  is  9.31xlOx  N/m2,  Lame  constant  p  is  1.03xl0x  N/m2,  fluid  density  p f  is  1000  kg/m3. 
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fluid  compressional  wavespeed  c  f  is  1500  m/s,  transducer  head  mass  M h  is  2.5  kg/m, 
transducer  tail  mass  MT  is  10.0  kg/m,  transducer  stiffness  Ks  is  lxlO6  N/nr,  transducer 
separation  distance  L  is  0. 1  m,  transducer  face  width  w  is  0.095  m,  transducer  stack  height  t  is 
0.01  m,  and  transducer  constant  g.33  is  0.025  (volts  m)/N.  In  figure  3,  the  solid  line  is  the  elastic 
plate  theory  developed  in  sections  2  through  4  and  corresponds  to  equations  ( 1 6)  through  (44); 
and  the  x  symbols  are  the  Bemoulli-Euler  plate  theory  and  correspond  to  equations  (45)  through 
(47).  The  elastic  plate  model  was  calculated  using  eleven  modes  (-5  <  n  <  5)  that  produced  a 
20-by  20-element  system  matrix  while  the  thin  plate  model  was  calculated  using  5 1  modes 
(-25  <n<  25).  Note  that  there  is  agreement  between  the  two  models  over  the  entire 
wavenumber  region. 
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Phase  Angle  (degrees)  Magnitude  (dB  ref  Volts /jiP a) 
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Wavenumber  (rad/m) 

Figure  3.  Transfer  Function  of  Transducer  Voltage  Divided  by  Input  Pressure 
Versus  Wavenumber  at  a  Frequency  of  100  Hz 
( _ )  Elasticity  Theory  and  (  x)  Thin  Plate  Theory 
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6.  A  NUMERICAL  EXAMPLE 


A  numerical  example  to  illustrate  the  dynamics  of  a  sonar  window  interacting  with  the 
Tonpilz  transducer  array  is  now  presented.  To  understand  fully  the  model  results,  it  is  necessary 
at  this  point  to  slightly  digress.  The  features  present  in  the  Tonpilz  transducer  output  are  best 
understood  if  the  dispersion  curve  for  the  system  is  studied.  The  dispersion  curve  is  found  by 
setting 


det  [A  -  F]  =  0 , 


(48) 


and  this  determines  the  location  in  the  wavenumber-frequency  plane  where  free  wave 
propagation  can  exist.  Each  of  these  free  waves  is  related  to  a  specific  dynamic  motion  of  the 
system.  Figure  4  is  the  dispersion  curve  for  an  example  with  the  following  system  parameters: 
window  thickness  h  is  0.1  m,  window  density  p  is  1200  kg/m3,  Lame  constant  A  is  1.43x10'' 
N/m2,  Lame  constant  p  is  3.57x1 08  N/m2,  fluid  density  p  f  is  1000  kg/m3,  fluid  compressional 
wavespeed  c  f  is  1500  m/s,  transducer  head  mass  Mh  is  1  kg/m,  transducer  tail  mass  Mj  is 
4  kg/m,  transducer  stiffness  Ks  is  1  x  1 07  N/m2,  transducer  separation  distance  L  is  0.1  m, 
transducer  face  width  w  is  0.095  m,  transducer  stack  height  /  is  0.0 1  m,  and  transducer  constant 
g33  is  0.025  (volts  m)/N.  The  noticeable  feature  from  figure  4  is  that  the  system  is  very  rich  with 
free  wave  propagation  above  2490  Hz.  This  frequency  value  corresponds  to  the  first 
antisymmetric  Lamb  wave  of  the  system.  Without  the  mass-loading,  this  would  occur  at 

f  =  —  =  2730  Hz  ;  (49) 

2  h 

however,  because  the  transducers  mass-load  the  plate,  this  location  is  shifted  downward  in 
frequency.  It  is  noted  that  the  multiple  free  waves  correspond  to  higher  order  plate  waves, 
fluid/structure  interaction  waves,  and  the  first  set  of  periodic  waves,  i.e.,  n  =  ±1.  These  periodic 
waves  are  related  to  the  spacing  of  the  individual  transducers  and  occur  in  wavenumber  at 
integer  multiples  of 


k 


2  n 
~L' 


(50) 
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Frequency  (Hz) 


Figure  4.  Dispersion  Curve  Versus  Arrival  Angle  and  Frequency 
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These  periodic  waves  are  frequently  called  Floquet  and/or  Bloch  ’  waves.  An 
additional  note  to  the  dispersion  curve  is  included.  Although  equation  (48)  may  predict  a 
specific  free  wave,  it  may  not  necessarily  propagate  in  the  structure.  (An  acoustic  load  may  not 
excite  a  specific  free  wave,  or  the  excitation  may  result  in  an  extremely  small  response.)  The 
propagation  of  a  wave  is  dependent  on  the  interaction  of  the  structural  load  with  the  physics  of 
the  free  wave. 

Figure  5  is  a  color  image  of  the  transfer  function  of  transducer  voltage  divided  by  input 
pressure  versus  arrival  angle  and  frequency  for  the  example  problem.  The  scale  is  shown  next  to 
the  image  as  a  colorbar  and  corresponds  to  dB  referenced  to  volts/pPa.  The  response  of  the 
system  shown  in  figure  5  is  extremely  complex,  however,  most  of  the  features  that  are  present 
are  nulls.  To  understand  the  relationship  of  these  nulls  with  respect  to  the  free  wave  propagation, 
the  dispersion  curve  is  plotted  over  the  transfer  function  in  figure  6.  The  color  map  scale  in 
figure  6  is  the  same  as  that  of  figure  5.  Nulls  frequently  exist  between  branches  of  the  dispersion 
curve  that  have  split  or  separated.  This  separation  is  due  to  the  fluid  loading  on  the  plate.  Figure 
7  is  four  cuts  of  figure  5  at  1500,  2610,  4030,  and  7320  Hz,  respectively,  and  is  shown  to 
illustrate  the  magnitude  (solid  line)  and  phase  angle  (dashed  line)  differential  versus  arrival  angle 
at  various  frequencies.  The  last  three  frequencies  are  shown  to  illustrate  the  nulling  effects  (or 
system  zero  dynamics)  that  are  present  at  higher  frequencies. 

Finally,  the  system  response  can  be  beamformed  using 


N 

B(k,a>)  =  ^ RVS{k,co )  exp[i(&  -ks)xn], 

n= 1 


(51) 


where  N  is  the  number  of  sensors  in  the  array,  ks  is  the  steered  wavenumber  (rad/m),  and  xn  is 
the  location  of  the  wth  sensor  (m).  For  the  following  analysis,  a  16-sensor  (element)  array  is  used. 
The  beam  patterns  are  displayed  as  polar  plots,  with  the  solid  lines  in  the  plots  corresponding  to 
array  response  based  on  the  theory  developed  previously  (equations  ( 1 )  through  (44)),  and  the 
dashed  lines  are  the  response  of  the  array  to  unity  input  at  all  wavenumbers.  Each  short  dashed 
concentric  half  circle  represents  10  dB  of  energy.  Because  these  beam  patterns  are  being 
displayed  as  polar  plots,  the  wavenumbers  in  equation  (51)  have  been  converted  into  arrival 
angles.  Figure  8  is  a  plot  of  the  beamformed  response  at  1500  Hz  with  a  steer  angle  of  0°. 

Figure  9  is  a  plot  of  the  beamformed  response  at  1500  Hz  with  a  steer  angle  of  30°.  It  is  noted 
that  in  these  cases,  the  sonar  window  improves  the  response  of  the  beamformer.  This  is  primarily 
due  to  the  receive  energy  drop  off  at  large  arrival  angles.  Figure  10  is  a  plot  of  the  beamformed 
response  at  26 1 0  Hz  with  a  steer  angle  of  0°.  Again,  the  sonar  window  dynamics  have  improved 
the  beam  pattern.  Figure  1 1  is  the  beamformed  response  at  2610  Hz  with  a  steer  angle  of  45.1°, 
an  angle  that  corresponds  to  a  null  in  the  system  response.  Steering  the  array  response  into  this 
null  location  has  significantly  degraded  the  response  of  the  beamformer.  Figure  12  is  a  plot  of  the 
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beamformed  response  at  4030  Hz  with  a  steer  angle  of  0°.  In  this  case,  the  sonar  window  slightly 
degrades  the  beamformer  response.  Figure  13  is  a  plot  of  the  beamformed  response  at  4030  Hz 
with  a  steer  angle  of  12.5°.  Steering  into  the  null  does  not  seriously  degrade  the  beamformer 
response.  Figure  14  is  a  plot  of  the  beamformed  response  at  7320  Hz  with  a  steer  angle  of  0°. 

The  sonar  window  significantly  improves  the  sidelobe  response  of  the  array  at  this  frequency  and 
steer  angle.  Figure  15  is  a  plot  of  the  beamformed  response  at  7320  Hz  with  a  steer  angle  of  9.8°. 
In  this  case,  the  beam  pattern  has  been  mildly  degraded. 
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Figure  5.  Transfer  Function  of  Transducer  Voltage  Divided  by 
Input  Pressure  Versus  Arrival  Angle  and  Frequency 
(Units  of  dB  re  volts/ /jP a) 
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Figure  6.  Dispersion  Curve  Plotted  Over  the  Transfer  Function 
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Figure  7.  Transfer  Function  of  Transducer  Voltage  Divided  by  Input  Pressure  Versus 
Arrival  Angle  at  (a)  1500  Hz,  (b)  2610  Hz,  (c)  4030  Hz,  and  (d)  7320  Hz 
( _ )  Magnitude  and  ( - )  Phase  Angle 
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Figure  8.  Beamformed  Response  at  1500  Hz  with  a  Steer  Angle  ofO° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  9.  Beamformed  Response  at  1500  Hz  with  a  Steer  Angle  of  30  ° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  10.  Beamformed  Response  at  2610  Hz  with  a  Steer  Angle  ofO° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  11.  Beamformed  Response  at  2610  Hz  with  a  Steer  Angle  of  45. 1  ° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  12.  Beamformed  Response  at  4030  Hz  with  a  Steer  Angle  ofO° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  13.  Beamformed  Response  at  4030  Hz  with  a  Steer  Angle  of  12.5  ° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  14.  Beamformed  Response  at  7320  Hz  with  a  Steer  Angie  of0° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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Figure  15.  Beamformed  Response  at  7320  Hz  with  a  Steer  Angle  of  9.8° 

( _ )  Tonpilz  Array  with  Window  and  ( _ )  Tonpilz  Array  without  Window 
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7.  CONCLUSIONS 


The  dynamics  of  a  fully  elastic  sonar  window  interacting  with  an  array  of  single  resonant 
Tonpilz  transducers  have  been  analytically  modeled.  This  model  shows  the  effects  of  acoustic 
energy  insonifying  the  array  at  all  frequencies,  rather  than  at  a  single  frequency.  This  is  useful  as 
sonar  systems  are  built  or  modified  for  broadband  processing.  The  model  predicts  where  the 
acoustic  cone  is  smooth  for  optimum  sonar  processing  and  where  the  effects  of  null  responses 
will  enter  into  the  processing.  The  effect  of  the  sonar  window  on  the  beam  pattern  has  been 
demonstrated. 
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APPENDIX 

MATRIX  AND  VECTOR  ENTRIES 


The  entries  of  the  matrixes  and  vectors  in  equation  (44)  are  listed  below.  Without  loss  of 
generality,  the  top  of  the  plate  is  defined  as  z  =  b  =  0.  For  the  [A("*(£  +  2m!  L )]  matrix,  the 


nonzero  entries  are 

2 

2  2  2  P  f 

a\  i  =  anX  2a  np  Akn  +  , 

Yn 

(A- 1 ) 

2 

1  9  9  CZ P  f 

a\ 2  =  a2  A  lain  Ak 2 

Yn 

(A-2) 

k„co2pf 

a\3  ~  2l*knPn  > 

Yn 

(A-3) 

knco2  Pf 

a\A  ~  2t^knPn  ’ 

Yn 

(A-4) 

a2\  =  — , 

(A-5) 

a22  =  2^nan  > 

(A-6) 

°23  =  Mifin  ~k„)  > 

(A-7) 

°24  =  MA7  “*«)  » 

(A-8) 

a3l  =(-a^A-2«^//-/W:^)exp(ia„a)  , 

(A-9) 

a32  =  (-a2A-2a2p-AJc2)exp(-ia„a)  , 

(A- 10) 

A- 


«33  =  -2 L&nPn  exp(i pnd)  , 


(A- 1 1 ) 


«34  =  2 fjknpn  exp(-i P„a)  , 

(A- 12) 

a4]  =  -2juknan  exp(i«„a)  , 

(A- 13) 

a42  =  2juknan  exp(-\ana)  , 

(A- 14) 

«43  =  /*(#?  -£/7)exp(i/?„a)  » 

(A- 15) 

and 

044  =  Mi  fin  -^)exp(-i/?„a)  , 

(A- 16) 

where 

kn=k+  l  , 

(A- 17) 

an  ~  ^  ~  ktl  , 

(A- 18) 

Pn=^~k2n  , 

(A- 19) 

and 

rn=J(«>/cf)2  -kl  =^k2f-kl  . 

(A-20) 

For  the  [F("*(£  +  2 nnl L)]  matrix,  the  nonzero  entries  are 

F  |  co ) 

h\  =  2 L  (ia„)exp(iawfl)  , 

(A-2 1 ) 

A-2 


hi  =-^^(-i«„)exp(-iawa)  ,  (A-22) 

hi  =  — (i^w)exp(i Pna)  ,  (A-23) 

/34  =  ~~~  (i^w ) exp(-i/7wa)  ,  (A-24) 

/4l  =  (i*„)exp(ia„a)  ,  (A-25) 

hi  =  Fjf^(^/7)exp(-iawa)  ,  (A-26) 

Lj 

hi  =  Fx-~- (~i/3n ) exp( \fln a)  ,  (A-27) 

Lj 

and 

/44  =  F-v^  (i/?„  )exp(-i pna)  .  (A-28) 

The  y  vector  entries  are 

y={.„  A'-'>  fi«-»  C<-‘>  O''1*  Jt0}  B<°> 

J  *  (A-29) 

C(°)  D(0)  ^(1)  g<i)  C(I>  0(1)  y 

The  p  vector  entries  are 

p  =  {-  2  Pj{(o)  0  0  0}T.  (A-30) 


A-3  (A-4  blank) 
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